(* Quadratura de f composta de ponto medio,usando particao uniforme em n subintervalos de [a,b];
precisao casas decimais: prec *)

f[x_]:=Sin[x];
a=0;
b=Pi;
Nmax=10;
prec=10;

Print["n valor quadratura"]; 
For[n=1,n≤Nmax,n++,
Print[n," ",
N[(b-a)/n Sum[f[a+(i+1/2)(b-a)/n],{i,0,n-1}],prec]
]
];